Fronts Propagation at the Onset of Plastic Yielding 
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The existence of a finite threshold, the yield stress, for the onset of plastic yielding is a universal 
feature of plasticity. This jamming-unjamming transition is naturally accounted for by the dynamics 
of a bistable internal state field. We show, within the athermal Shear Transformation Zones (STZ) 
theory of amorphous plasticity, that the transition is accompanied by the propagation of plastic 
fronts. We further show that the mean-field theory cannot select the velocity of these fronts and go 
beyond the mean-field description to include fluctuations and correlations effects, resulting in new 
nonlocal terms in the equations. Finally, we demonstrate that the new terms, with an associated 
intrinsic lengthscale, provide a velocity selection mechanism for the plastic fronts. 
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I. INTRODUCTION 

A complete understanding of the dynamics of plastic 
deformation in low temperature, or athermal, amorphous 
systems remains a major theoretical challenge in statisti- 
cal and condensed matter physics. The response of these 
systems to the application of external driving forces ex- 
hibits some universal features like a transition from jam- 
ming to flow at the yield stress, orientational memory 
effects of recent deformation and strain localization 
In this work we focus on the spatiotemporal patterns as- 
sociated with the transition from jamming to flow at the 
onset of plastic yielding, in the framework of the recently 
developed athermal Shear Transformation Zones (STZ) 
theory of amorphous systems 0, E|. This transition is 
controlled by the applied stress. For stresses below a 
material dependent threshold, the yield stress, these sys- 
tems exhibit elastic deformation as well as transient plas- 
tic flow that vanishes at a finite time. The resulting state 
is "jammed" in the sense that it carries a finite stress 
without flowing. On the other hand, for stresses that 
exceed the yield stress, these systems exhibit persistent 
plastic flow, implying some remarkable "unjamming" (or 
yielding) transition at the yield stress. 

The main difficulty in developing a theory of amor- 
phous plasticity, that among other things should predict 
plastic yielding at a finite yield stress, is that the iden- 
tity and characteristics of the microstructural objects 
that "carry" plastic deformation in amorphous (disor- 
dered) systems remain quite elusive. However, a growing 
body of experimental and simulational evidence seems to 
point towards a unifying picture in which low tempera- 
ture amorphous plasticity involves stress driven config- 
urational rearrangements of localized regions composed 
of a small collection of the relevant elementary entities 
(e.g. particles, molecules, grains, colloids, bubbles) 
The existence of such a unifying picture, which is inde- 
pendent of the elementary entities and their microscopic 
interactions, is fully consistent with the existence of uni- 
versal features mentioned above. 

In line with this accumulating evidence, the athermal 
Shear- Transformation- Zones (STZ) theory of amorphous 



plasticity views STZ as localized groups of particles that 
are more susceptible to shearing transformation under 
stress than their surroundings. Upon surpassing a lo- 
cal threshold an STZ can undergo a finite irreversible 
shear in a given direction. Once transformed, due to 
a local redistribution of stresses, an STZ resists further 
deformation in that direction, but is particulary sensi- 
tive to reverse shear deformation if the local applied 
stress changes sign. Therefore, an STZ is conceived as 
a two-state system that can transform between its in- 
ternal states depending on the magnitude and direction 
of the local stress. In addition, the stress redistribution 
that accompanies an STZ transition can induce the cre- 
ation and annihilation of other STZ at a rate proportional 
to the local energy dissipation (recall that thermal fluc- 
tuations are assumed to be absent or negligible). The 
minimal two-state assumption, along with the deforma- 
tion driven creation and annihilation of STZ, provides a 
mechanism for retaining and losing orientational memory 
of previous deformation and is sufficient to capture the 
transition from jamming to flow at a finite yield stress. 

We choose to study the spatiotemporal characteristics 
of the unjamming transition using STZ theory for two 
major reasons. First, the existence of the transition is 
described in this theory, as will be explained in detail be- 
low, by an exchange of dynamic stability in the equations 
of motion for an internal state field. We find such a de- 
scription superior to theoretical approaches that simply 
postulate the existence of a threshold stress above which 
plastic flow persists More generally, we find an ap- 
proach that aims at identifying proper order parameters, 
or internal state fields, and writing equations of motion 
for them based on microscopic insights as well as symme- 
tries and conservation laws very appealing 0, Q • Second, 
recent work has shown that the STZ equations capture a 



variety of phenomena observed in computer simulations 
and laboratory experiments @, @, & ©, E3, ED, GJ, El El • 
surposes here are the results re- 



Most relevant for our 

ported in [l3|, E3 . In [13| , the two-state nature of STZ 
has been observed directly in bubble-raft experiments 
and two-dimensional foam simulations. It was found that 
groups of bubbles can make transitions, the so-called Tl 
events that are realizations of STZ transitions in foam, 
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between two states if the direction of the applied force is 
reversed shortly after an event occurs. However, if the 
force is not reversed until after other events have occurred 
nearby, then the orientational memory is lost, implying 
that an STZ is annihilated. In [lj] the STZ equations 
were shown to agree quantitatively with transient shear 
reversal behavior of a granular flow in a Taylor-Couette 
cell. This result provides the first direct support to an 
equation that will be central in the analysis to follow. 

We proceed by introducing the relevant mean-field STZ 
equations in Sect. [Til where we also discuss the existence 
of the transition from jamming to flow at a finite yield 
stress and the physics behind it. In Sect. |TTT]we focus 
on the spatiotemporal characteristics of the unjamming 
transition. We show that there exist propagating front 
solutions if the equations are modified to include effects 
beyond mean-field as well as an intrinsic lengthscale. Sec- 
tion IIVI offers some concluding remarks. 



II. ELEMENTS OF THE ATHERMAL SHEAR 
TRANSFORMATION ZONES (STZ) THEORY 

Consider a two-dimensional athermal amorphous sys- 
tem under the application of a pure shear stress s. The 
two states of an STZ are oriented along the principal axes 
x (the direction of the force) and y (the perpendicular di- 
rection) , and denoted by "±" . The number densities of 
STZ in the "±" states are denoted by n±. The configu- 
rational disorder of the system is characterized by an ef- 
fective disorder temperature \ that is assumed to govern 
the total density of STZ, with n^e -1 /* being the steady 
state density of STZ under persistent deformation, x was 
shown very recently to play an important role in the de- 
formation dynamics of amorphous solids @, H, H, QU [H| ; 
however, its dynamics in the present context are not of 
prime interest and will not be discussed. The internal 
state fields n± and x> m addition to the applied stress 
s, determine the macroscopic plastic rate of deformation 
D pl according to [U 



T D pl = \(R{s)n--R{-s)n^j , 
tq n± = R(±s)n T — -R(T S ) n± + T 
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Here A is a material-specific parameter with the dimen- 
sions of area, tq is the basic time scale and R(±s)/to 
are the rates for forward and backward transitions re- 
spectively. The athermal condition implies that R(s) 
vanishes for s < 0, i.e. there are no transitions in the 
direction opposite to the direction of the applied force 
1- 

Eq. (UJ) simply states that the macroscopic plastic rate 
of deformation D pl results from STZ transitions between 
their two internal states, depending on the driving force 
through R(s) and on the population of STZ in each of 
the two states. Eqs. © are Master equations for the 
population densities n± themselves. The first two terms 



on the right hand side describes the number conserving 
process in which STZ change their state from "+" to "-" 
or vice versa. The third term on the right hand side is 
a coupling term that accounts for interactions between 
STZ. As mention above, the stress redistribution that 
accompanies an STZ transition can induce the creation 
and annihilation of other STZ . This is a crucial aspect 
of the theory. The rate of STZ creation T, which is a 
positive definite scalar, is assumed to be proportional 
to the rate of plastic dissipation 2sD pl 0, H| . For low 
temperature we obtain [3, |3j 
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where s y is a material-specific parameter of stress dimen- 
sion, to be shown below to be equal the yield stress. s y 
has an interesting physical meaning: it quantifies how 
much of the dissipated energy is invested in creating new 
STZ. When it is large, only a small amount of the dissi- 
pated energy is invested in creating new STZ (available 
for irreversible transitions) and vice versa. This inter- 
pretation immediately suggests that s y is related to the 
yield stress; see [l5[ for a discussion of this issue from an 
entropic point of view. 

We proceed by setting the total STZ density to its 
steady state value, n++n_ =nooe _1 ^ x and by defining 
an orientational order parameter m (which is generally a 
tensor, see Q) as 



(4) 
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This internal state field represents the bias of the STZ 
populations between the ± states. With these choices we 
can rewrite Eqs. ©-(HI) as 



T Q D pl = e e- 1/x C(s) 
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_ R(s) + R(-s) 
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The theoretical framework in which the amorphous sys- 
tem is characterized by coarse-grained internal state 
fields, in addition to stress and strain, is a major point of 
divergence compared to other approaches Q and plays a 
crucial role in the discussion below, where the field m is 
the central object. 

In [3], Eqs. (J5|)-([6|) were directly applied to predict 
transient granular flow in experiments of shear reversal 
in a Taylor-Couette cell. The theory was shown to agree 
well with the experimental data, encouraging us to fur- 
ther study it in the present context. The first issue to be 
discussed is the stable steady state solutions of Eq. © . 
It has two stationary solutions: 
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1. A jammed state with m = ±1 (depending on the 
sign of s), where all the STZ arc cither in the "+" 
or the "-" state respectively, and D pl = in Eq. 
©> 

2. A flowing state with m = s y /s and D pl ^0 in Eq. 
©■ 

These two solutions coincide at s/s y = ±l. It is straight- 
forward to show that the jammed state is dynamically 
stable for \s\ < s y , while the flowing state is dynami- 
cally stable for \s\ >s v . An exchange of stability occurs 
at s = Sy, which indeed shows explicitly that s y is the 
yield stress. The fact that s y was introduced in Eq. © 
as a proportionality constant between the rate of STZ 
creation T and the rate of energy dissipation 2sD pl is 
yet another point of divergence with respect to standard 
approaches. In STZ theory the yield stress s y , as was ex- 
plained before, is intimately related to the ability of the 
material to create new STZ as a result of plastic dissipa- 
tion, i.e. STZ transitions in nearby locations [TEl ]. In this 
interpretation, a material can yield, i.e. support steady 
plastic flow, if enough STZ per unit time are created due 
to other STZ transitions. If this rate is not sufficiently 
large, then existing STZ are exhausted in transitions at 
a given direction and not enough new STZ are available 
for additional transitions that are needed for attaining a 
steady state. For a given applied stress, whether or not T 
is large enough, is determined by s y . This is the essence 
of yielding in STZ theory. 

As a result, the onset of plastic yielding at a finite 
threshold is a dynamic phenomenon manifested as an ex- 
change of stability of the bistable field m @ • These re- 
sults have similar implications for strain-rate 7 controlled 
experiments (complementary to the stress-controlled ex- 
periments discussed up to now) where the existence of 
a finite yield stress means that a steady state with a fi- 
nite stress is obtained in the limit 7 — > 0. The existence 
of a dynamic threshold is a signature of an underlying 
bistability (or multistability), which is a natural way to 
get a finite dissipation at a vanishing strain-rate. Note 
that all these conclusions are entirely independent of the 
material function C(s). We thus propose that the bista- 
bility of the orientational order parameter m is a crucial 
feature of the theory of amorphous plasticity. 

III. THE TRANSITION FROM JAMMING TO 
FLOW 

We now proceed to analyze the spatiotemporal char- 
acteristics of this jamming- unjamming transition. We 



substitute Eq. ([5]) into Eq. ([6]) to obtain 

T m = 2C{s){l-m)[l-—J , (8) 

where we have assumed s > 0. Consider then an exper- 
iment in which the stress is ramped to a constant value 
smaller than the yield stress (i.e. s < s y ). The system 
is assumed to be initially undeformed and isotropic such 
that the orientational order parameter satisfies m(t = 0) = 
0. Eqs. © and ([H]) then predict that sub-yield plastic de- 
formation occurs as m relaxes to the jamming fixed-point 
m = 1 on a typical time scale Tj am (s) 



during which D pl — > 0. The jamming time scale Tj am (s) 
diverges as s—+s~ , i.e. approaching the yield stress from 
below, and as s — > + . The latter divergence is a result 
of the athermal condition, leading to C(s — > 0) — > 0|. 
Therefore, Eq. ([9]) predicts that jamming is obtained 
by progressively slower creep-like sub-yield deformation 
as s approaches either or 1. Moreover, it provides a 
way to measure the phenomenological function C(s) in 
the range where the jamming time is experimentally (or 
simulationally) accessible. 

Suppose now that s/s y is 0(1), but not very close to 
unity such that the system becomes jammed on a realistic 
time scale; then suppose that after jamming the applied 
stress is ramped again to a value above the yield stress, 
i.e. s > s y . The system is now in the jammed state, 
m = 1, that is unstable for s > s y . In these situations 
we generally expect fluctuations to propagate as fronts, 
converting an unstable state into a stable one. Therefore, 
we are looking for translational invariant solutions of the 
form m(x,t) = m[x — ut)\ here we assume that the x- 
dimension is much larger than the y-dimension and that 
the front propagates from left to right. Our task now 
is to solve Eq. |(5J| for this ansatz, with the boundary 
conditions 

g 

m(x—>—oc)—t — , mix — > +oo)^l . (10) 
s 

The first boundary condition corresponds to the flowing 
state left behind the front, while the second one corre- 
sponds to the jammed state ahead of it. A solution can 
be readily obtained, yielding 
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m(x, t) = 



exp 



-2C(s)(s/s y - l)(x - ut) 



T U 



— 1 ) exp 



-2C(s)(s/s y - ut) 
t u 



(11) 



r 



This seems to be a propagating front solution in which a 
plastically deforming region, D pl ^ 0, invades a jammed 
region, D pl = 0. 

However, Eq. (jlip is not a unique solution of the prob- 
lem since it is valid for any velocity u. More importantly, 
a velocity u cannot be selected by Eq. ([8]) in principle 
since in the absence of an intrinsic lengthscale a veloc- 
ity cannot be dimensionality constructed. Therefore, if 
we aim at describing the spatiotemporal patterns that 
accompany the onset of plastic flow, we cannot avoid ad- 
dressing the fundamental problem of missing a length- 
scale in our theory. In fact, two lengthscales are al- 
ready implied within our theoretical framework; a fi- 
nite STZ density implies a typical distance between STZ, 
which can be thought of as a microstructural correlation 
length. This lengthscale might be relevant as different 
STZ can interact via the stress and displacement fluc- 
tuations generated when an STZ transforms between its 
internal states. An evidence for the existence of such a 
lengthscale was given, for example, in [H, [13; it can be 
identified with the typical size of areas of quasi-crystal- 
like short-range order of [l6| or with the typical size of 
regions free of liquid- like defects of [l7| • 

The range in which these interactions can affect STZ 
transitions also implies a lengthscale, possibly larger than 
the previous one, that is determined by a combined ef- 
fect of the magnitude of the stress fluctuation, its range 
and the distribution of STZ transition thresholds. This 



lengthscale characterizes the scale in which the probabil- 
ity of finding a sufficiently large stress perturbation (that 
can overcome the local transition threshold) is not expo- 
nentially small. Therefore, we should consider nonlocal 
stress fluctuations that can induce STZ transitions at dif- 
ferent locations where STZ already exist and are close 
to their transition threshold. This effect was stressed 
by several authors [la, [13, |20| and was assumed to be 
the origin of the cascades of rearrangements observed in 
athermal quasi-static simulations [2l[ . With this physical 
intuition in mind, we should extend the homogeneous and 
local mean-field theory to include some inhomogeneous 
and nonlocal effects. 

In summary, our theory should incorporate nonlocal 
terms that account for the following physical effects: (i) 
The stress and displacement fluctuations that are gener- 
ated by STZ transitions, (ii) The joint spatial probability 
distribution function of STZ that includes a correlation 
length, (iii) The STZ transition thresholds distribution. 
Therefore, we add to the left hand side of Eqs. (HI)-© 
terms that are weighted integrals over existing terms that 
include the transition rates i?(±s); these new terms rep- 
resent that idea that STZ transitions from "+" states to 
"-" states (or vice versa) at a given location result in a 
properly weighted change in the rate of transitions at dif- 
ferent locations. Rearranging the modified equations we 
arrive at the following equation for m(x,t) 



T$m(x, t) = 2 C(s) [1 — mix, t)] 



s m(x, t) 



2 C(s) 



1 - 



s m(x, t) 



K(x-x')[l-m(x',t)]dx'. (12) 



r 



The phcnomcnological kernel K(x) in Eq. lj!2|L whose 
dimension is {length)" 1 , represents the physics of STZ 
correlations and interactions discussed above. The fact 
that it is a one-dimensional function of x alone already 
incorporates several approximations. First, there is some 
evidence that STZ transitions generate quadrupolar elas- 
tic fields [HI] . Here this anisotropic structure is neglected 
due to the assumption that the y-dimension is much 
smaller than the x-dimension. Second, the kernel K 
might include some time dependence that accounts, for 
example, for the wave nature of the stress and displace- 
ment fluctuations. We assume that this time-dependence 



can be integrated out without affecting the basic physics 
we are trying to describe. In choosing a specific func- 
tional form for the kernel K[x) we do not want to make 
any definite claim as to whether the incorporated length- 
scale is a short-distance correlation length or a longer- 
distance cutoff scale discussed above, but to stress that 
our theoretical considerations entails the existence of a 
lengthscale. Both physical possibilities share the feature 
that the immobile STZ interact via nonlocal stress flue- 
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tuations which can be schematically described by 



(13) 



where £ is the characteristic lengthscale and a is the am- 
plitude. Note that we do not consider here very long- 
range kernels, i.e. with power-law tails, as we believe 
that the STZ transition thresholds distribution limits the 
range in which stress fluctuations can induce STZ tran- 
sitions, as explained above. 

Eq. (fT2)l . with the kernel of Eq. (fT3|) , poses a non- 
trivial front selection problem. As in many other front 
selection problems [22[, there might exists a family of 
propagating front solutions from which a unique solution 
is selected dynamically. A complete analysis of Eq. (TT2]) 
is well beyond the scope of the present work. However, 
we aim here at demonstrating that the proposed new 
terms indeed lead to front selection. To that aim, we first 
consider the possibility that the nonlocal contribution is a 
small correction to the mean-field equation and treat Eq. 
(I12p perturbatively with a being the small parameter in 
the expansion. We thus write m — m^+amW+0(a 2 ), 
where to/ ) is given in Eq. (jlip . The first order equation 
in a selects the velocity, from which we estimate 
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FIG. 1: (Color online) The propagating m front solution of 
Eq. (|12|l with a — 0.01 and s = 1.1 (solid line) compared 
to the analytic prediction given by Eq. (|lip (dashed line), 
shifted by a few time units for clarity. (Inset) The velocity u 
(in units of 21/ to) as a function of s, both for the numerics 
(open circles) and the prediction of Eq. (|14|l (solid line). 



IV. CONCLUDING REMARKS 



It is important to note that the nonlocal term is a singular 
perturbation in the sense that there is no solution at all 
in the limit 

As we cannot offer a similar analysis for the general 
case where a is not necessarily small, we studied Eq. 
(|12p numerically by choosing s> s y and introducing the 
homogeneous jammed state, m = l, with various localized 
perturbations. In the numerical calculations we choose, 
for simplicity, C(s)=H(s — s y )(s — s y ) for s > 0. Here H(-) 
is the Heaviside unit step function. We have found that 
all perturbations converge, for sufficiently large times, to 
a unique front profile with a unique velocity u. This result 
shows that indeed the new nonlocal term in Eq. (jX3J) 
provides a front velocity and profile selection mechanism 
for the problem at hand. In Fig. [JJwe show the numerical 
front profile for a small amplitude a and the analytic 
prediction of Eq. (fTTj) (shifted by a few time units for 
clarity). In the inset we compare the numerical results 
for the velocity u as a function of s to the prediction 
of Eq. (dU). We observe that in this regime, i.e. for 
small a, the profile is, to a very good approximation, 
given by Eq. (fTTj) with the selected velocity satisfying 
Eq. HU). Upon increasing the amplitude a, the nonlocal 
term becomes more dominant and the solution of Eq. 
([Til is no longer accurate. We thus conclude that in 
the presence of the nonlocal interaction term, with the 
associated lengthscale £, the proposed theory predicts the 
existence of propagating front solutions that remain to be 
observed in computer simulations and experiments. 



In this paper we asked what are the spatiotempo- 
ral characteristics of the transition between a jammed 
state found at applied stresses below the yield stress 
s y and a homogeneously plastic flow state found at ap- 
plied stresses above the yield stress s y . In our opinion, 
this is a fundamental question in the field of plasticity 
of amorphous systems. We addressed this problem in 
the framework of the recently developed athcrmal Shear- 
Transformation-Zoncs (STZ) theory, which in our opin- 
ion offers a promising route for developing a predictive 
theory of amorphous plasticity. This theory describes 
the phenomenon of plastic yielding at a finite thresh- 
old in terms of an exchange of dynamic stability in the 
equation of motion for an orientational order parameter 
to. Our main result is that this theory predicts the exis- 
tence of propagating front solutions that accompany the 
yielding transition provided that nonlocal effects are in- 
troduced. The new nonlocal terms, that account for non- 
local STZ interactions, incorporate a lengthscale that is 
missing in the original theory. Our predictions for the ex- 
istence of plastic fronts at the yielding transition should 
be tested in computer simulations and experiments. This 
may provide further support for the analytic structure of 
the developing STZ theory, substantiate the existence of 
the bistable orientational order parameter to and shed 
some light on the missing intrinsic lengthscale. As the 
existence of a yield stress is a fundamental property of 
materials, elucidating the spatiotemporal nature of the 
yielding transition is of major importance. 
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